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Because of a formal equivalence with the partition function of an Ising chain, the semiclassical 
traces of the quantum baker map can be calculated using the transfer-matrix method. We analyze 
the transfer matrices associated with the baker map and the symmetry-reflected baker map (the 
latter happens to be unitary but the former is not). In both cases simple quantum-circuit repre- 
sentations are obtained, which exhibit the typical structure of qubit quantum bakers. In the case 
of the baker map it is shown that nonunitarity is restricted to a one-qubit operator (close to a 
Hadamard gate for some parameter values). In a suitable continuum limit we recover the already 
known infinite-dimensional transfer-operator. We devise truncation schemes allowing the calcula- 
tion of long-time traces in regimes where the direct summation of Gutzwiller's formula is impossible. 
Some aspects of the long-time divergence of the semiclassical traces are also discussed. 



PACS numbers: 05.45.Mt 



I. INTRODUCTION 



The quantum baker's map [l| is a very useful test- 
bench for investigating quantum-classical correspondence 
issues in a variety of settings. Conceived as a model for 
studying the semiclassical limit of closed chaotic systems, 
it was first applied to analyze the random-matrix conjec- 
ture, the scarring phenomenon, Gutwiller's trace formula 
and the long-time validity of semiclassical approxima- 
tions (see, e.g., [243 )■ the last years "open" quan- 
tum bakers were constructed and employed for studying 
semiclassical aspects of the scattering problem [sj, e.g., 
the fractal Weyl law for the distribution of resonances 
i- 

The quantum baker also appeared in a variety of prob- 
lems of Quantum Information, Quantum Computation 
and Quantum Open Systems. It was noted that the 
quantum baker could be efficiently realized in terms of 
quantum gates [13] ■ A three qubit Nuclear Magnetic 
Resonance experiment was proposed 11 ill and then im- 
plemented (with some simplifications) |12l |. 

On the theoretical side, Schack and Caves 13 [showed 
that the Balazs-Voros-Saraceno quantum baker % [H can 
be seen as a shift on a string of quantum bits (in full 
analogy with the classical case), thus taking an impor- 
tant step towards generalizing the method of symbolic 
dynamics to the quantum case. At the same time, their 
circuit representation led naturally to a family of alterna- 
tive quantizations. This family of bakers was the subject 
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of several studies [l3 - [T6j . The ability of the baker fam- 
ily to generate entanglement was first studied by Scott 
and Caves 0| (see also [l3|)- Decoherent variants of the 
baker map were constructed by including mechanisms of 
dissipation and/or diffusion [l8l - l2(| . 

From the point of view of the structure of the quantum 
baker map, important results were recently obtained by 
Ermann and Saraceno 2l|, who, building upon previous 
work by Lakshminarayan and Meenakshisundaram [23|. 
generalized the Schack-Caves family, and showed that all 
quantum bakers are perturbations of a common kernel 
(the "essential" baker). 

The present paper focuses on an almost unexplored 
aspect of the semiclassical theory of the quantum baker 
map. Gutzwiller's approximate formula for the traces 
of the baker map is formally equivalent to the partition 
function of a finite Ising chain (with exponentially de- 
caying interactions and imaginary temperature). Thus, 
it can be evaluated with the standard transfer-matrix 
method. We have studied extensively the transfer matri- 
ces associated with the semiclassical traces of the baker 
map and the symmetry-reflected baker map. Our most 
remarkable finding is the existence of a qubit structure 
hidden in the semiclassical traces: the transfer matri- 
ces admit a quantum-circuit representation that is very 
similar to that found by Schack-Caves for the quantum 
baker. 

We know of two studies of the baker map which ap- 
plied the transfer matrix method (to the semiclassical 
long-time propagation of wave-packets Ml , to the analy- 
sis of periodic-orbit action correlations [23|)- Diff^erent in 
spirit, the present work aims at analyzing the transfer- 
matrix method in itself. 

The paper is organized as follows. First we discuss 
the general structure of the transfer matrices for the 



2 



baker map, exhibiting their quantum-circuit represen- 
tation. Though the transfer matrices are not unitary, 
nonunitarity is restricted to a one-qubit "gate" . Except 
for this fact, the transfer matrices resemble qubit quan- 
tum bakers (Sec. Ill|. 

Section Hill is devoted to spectral properties. We show 
that all eigenvalues lie inside the unit circle in the com- 
plex plane, or very close to it. The number of eigenval- 
ues close to the unit circle coincides approximately with 
the quantum dimension (for a suitable parameter range). 
This situation is very similar to that found in the infinite- 
dimensional transfer-operator method |6|, l24l ) . Indeed we 
show that the transfer operator arises as a suitable con- 
tinuum limit of the transfer matrix. 

In Sec. II VI we exhibit truncation schemes which permit 
the calculation of long-time traces. 

Some aspects of the long-time divergence of the semi- 
classical traces are assessed in Sec. El 

The baker map possesses a spatial symmetry. If one 
is interested in separating the spectrum/traces into sym- 
metry classes, then one must consider also the symmetry- 
related transfer matrices. Though similar to the matrices 
of Seciniin many respects, these matrices happen to be 
exactly unitary. They are studied in Sec. IVII 



Concluding remarks are presented in Sec. IVII 



II. BAKER MAP: TRANSFER MATRIX 

In order to make the paper self-contained we start 
by summarizing some basic information about the baker 
map. 

The dynamics of the classical baker map is schemat- 
ically depicted in Fig. [H It acts on the unit square 
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FIG. 1: The classical baker map. The dynamics is controlled 
by the fixed points (0,0) (for q < 0.5) and (1,1) [q > 0.5). 
The action of the map is to compress/stretch the coordinates 
p/q relative to the fixed points by a factor of two. Some 
examples are depicted. 

piecewise- linearly. Points with q < 0.5 are governed by 



the hyperbolic point at the origin. If g > 0.5 the fixed 
point at (1,1) rules. In any case the coordinates p/q rel- 
ative to the fixed point are compressed/stretched by a 
factor of two: 



q' = 2q~e, 



(1) 
(2) 



where e = [2g], the integer part of 2q. 

The quantum analogue of the classical mapping is a 
unitary operator acting on an even-dimensional Hilbert 
space Xi^M- I'^ the g-representation its matrix reads Q : 



B = G 



N 



G 



N/2 







G'Ar/2 



(3) 



where Gjv is the A^-dimensional antiperiodic Fourier ma- 
trix: 



(G 



N 



-2TTi{n+l/2)(m+l/2)/N 



N 



(4) 



with < m, n < — 1. For this abstract system the 
Planck constant coincides with the inverse of the dimen- 
sion, i.e., fi = l/(27r7V). 

Our main concern are the traces of the quantum baker, 
tr B^, for L = 1, 2, . . .. In the semiclassical regime of N 
large enough (for a given L) the traces can be approxi- 
mated by the Gutzwiller formula @ 



2^/2 ^ 



,2TTiNS{e) 



(5) 



where the sum runs over all periodic trajectories of length 
L of the classical map, indexed by the binary column 
vectors e = (eq, ei, . . . , el-i)*, with = 0,1. The cor- 
responding actions S{e) are quadratic functions of the 
binary symbols. 



5(e) e'Ae 



(6) 



The "coupling" matrix A is suitably expressed in terms 
of the matrix 5' of a cyclic shift, S ■ (oq, oi, . . .ul^iY — 
(ai, . . .,aL-i,aoY Hi]: 



■)L-1 



A 



L-1 
i=0 



(7) 



The trace formula for the baker map ^ can be derived 
as follows. Write tr i?^ as a sum of products of Fourier 
matrix elements using Eq. ([3]) . Approximate sums by in- 
tegrals. Extend the limits of integration (which in prin- 
ciple are finite) to ±oo. The remaining integrals, being 
Gaussian, can be done exactly, giving the Gutzwiller sum 
i- 

Equation ([S]) is a particular case of the general 
Gutzwiller trace formula for systems with a chaotic clas- 
sical limit. For such systems the trace formula lays a 
bridge between the quantum spectrum and the set of 
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classical periodic orbits. It constitutes the core of all 
semiclassical schemes for relating energy levels to classi- 
cal information, from the pioneer attempts [H, [27l[, t o 
the most recent, highly sophisticated developments p8| . 

The baker is very special in that all its periodic orbits, 
and properties thereof, are known analytically. So, in 
principle, the periodic-orbit sum ([5]) can be calculated 
for any L. However, because of the exponential growth of 
the number of periodic orbits with period, the brute-force 
summation is restricted to, say, L < 30. Remarkably, 
the method developed by Dittes, Doron and Smilansky 
does not suffer from such a limitation. These authors 
defined an infinite-dimensional integral operator, 



S Q- 



< q, < 1, whose traces give exactly the periodic-orbit 
sum of Eq. They showed that tr can be efficiently 
calculated from the eigenvalues of the matrix of W in the 
Fourier basis, after suitable truncation [gI. [23|. 

An alternative to the infinite-dimensional operator 
method relies on the formal identification of the periodic- 
orbit sum ([5]) with the partition function of an Ising 
chain, for a purely imaginary temperature. The cycli- 
cal nature of the coupling matrix A says that this is a 
circular chain, consisting of L spins. All spins interact 
among themselves, but, as interactions decay exponen- 
tially with distance, one may expect some computational 
benefit (without much loss of accuracy) by truncating the 
interaction to some given number r of closest neighbors, 
i.e., 



A: 



2^-1 ^ 2^ 

i=0 



(9) 



The exact periodic-orbit sum corresponds to setting r = 
L — 1 (no truncation at all), but, in principle, one can 
consider any truncation, even to first neighbors (r ~ 1). 

Once one has established the equivalence between the 
baker periodic-orbit sum and the Ising partition func- 
tion, the transfer-matrix method can be invoked. The 
first- neighbor case is explained in textbooks [2^. It 
seems that the many-neighbor case has not been explic- 
itly worked out in the literature, but some discussions 
exist |27|, l30[ . Anyway, even if not trivial, the general- 
ization can be carried out following the spirit of the one- 
neighbor case and, in the case of the baker-Ising, leads 
to 



2L/2 



g27ri7V5(e) 



2^-1 



trM^ 



the explicit expression for M being [25 
/ 1 1 



V2 








ao ai 
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' -J 



The elements ak are given by 



Uk = exp 



with < fc < 2'' - 1, and 



1 + — 



TrN2^ 
2^~I 



(11) 



(12) 



(13) 



Note that M is a complex matrix depending on the three 
parameters {N,L,r). Its dimension is 2*" x 2'', with r < 
L—1. When r — L—1 one recovers the exact semiclassical 
traces, i.e., Eq. (jlOp becomes an equality. 

A prefactor 1 / \/2 has been absorbed into the definition 
of M because in this way most of the spectrum of M lies 
close to the unit circle or inside it (see Sec. IIIII below) . 
So, the transfer matrix M becomes qualitatively similar 
to the semiclassical transfer operator of Refs. 0, 
come back to this point later). 

The transfer matrix approach transforms the calcula- 
tion of the Gutzwiller sum into the problem of obtain- 
ing the trace of the L-th power of the finite matrix M. 
From a numerical point of view this is advantageous only 
if truncation of AI is admissible. We defer the analy- 
sis of this question until Sec. [V] Now we concentrate on 
the transfer matrices themselves which, as we shall show, 
possess very interesting properties. 

The transfer matrix (jlip exhibits the structure of a ten- 
sor product which is conveniently displayed by switching 
to a qubit representation. This consists in identifying the 
"spatial" degree of freedom q with the tensor product of 
r two-level systems: 



\k) 



(14) 



Here k indexes the states of the g-basis and e.^ = 0, 1 is a 
label for the qubit basis states; they are related through 
the binary expansion 



k = eo + 2ei 



(15) 



In the qubit picture the matrix M can be broken up 
into the elementary gates shown in the quantum circuit 
of Fig. m The circuit acts on r qubits, starting with a 
downwards cyclic shift of all qubits, i.e.. 



(10) 5^|eo)®|ei)<^---«)| 



>(8)|eo)«'|ei)«'- • 



-2) ■ 
(16) 
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III. SPECTRAL PROPERTIES 



a/ 
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FIG. 2: Quantum circuit associated with the semiclassical 
trace of the L-th power of the baker map. The parameter r < 
L — 1 controls the level of truncation {r = L — 1 corresponds 
to no truncation at all). From left to right, the circuit is 
composed of a downwards shift of all qubits (S"''), a nonunitary 
one-qubit gate H, and a sequence of symmetric phase gates. 



A nonunitary gate acting on the first qubit follows, which 
bears some resemblance with a Hadamard gate. 



72 V 1 ^ 



with 



exp 



(17) 



(18) 



It must be noted that, as the dimension N is an even 
number, we always have x ^ —1, and H is never unitary. 
After the H gate one finds a single qubit phase gate. 



Poo (a) 



1 

e^" 



(19) 



Finally one has a sequence of two-qubit phase gates, act- 
ing symmetrically between the first and the fc-th qubit, 



10 
10 
10 



(20) 







with exponentially decreasing phases, (3 — 
a/2,a/4:,...,a/2\ 

Writing M as a circuit has helped us in identifying 
its basic constituents. In particular we recognize the 
one-qubit gate H as responsible for the deviation from 
unitarity. If we substitute H hy a, Hadamard gate, 
i.e., set 2: — —1 in pT|) . then the circuit of Fig. [2] 
acquires the typical structure of the baker family (21)]: 
an "essential-baker" block (formed by a shift and a 
single-qubit Fourier transform) followed by a "diffraction 
kernel" (phase gates). 



Figure [3] displays a typical transfer-matrix spectrum in 
the case that the dimension of the transfer matrix is much 
larger than the quantum dimension, i.e., 2^~^ ^ N. The 
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FIG. 3: Comparison of spectra. Open circles correspond to 
the quantum baker map (A^ = 46) and dots to the semiclas- 
sical transfer matrix for L = 11. 

spectrum can be coarsely divided into three parts. Ap- 
proximately N eigenvalues are localized close to the unit 
circle. Most of the remaining ones are concentrated in 
a disk of smaller radius, and there are some transitional 
eigenvalues spiraling out from the inner disk to the unit 
circle. The eigenvalues lying close to the unit circle can 
be put in almost one-to-one correspondence with the ex- 
act eigenvalues of the quantum baker ([3]). 

Increasing L while keeping TV fixed does not signifi- 
cantly change the almost unitary part of spectrum, but 
populates the region of small moduli (see Fig. 0]). This 
fact, when combined with the existence of eigenvalues 
with moduli larger than one, leads to the conclusion that 
Gutzwiller's traces for the baker are divergent in the limit 
L — > oo (fixed TV). This result has been confirmed previ- 
ously using the transfer-operator method @, [23 ? ] . 

Indeed, the overall features of the transfer- matrix spec- 
trum described above can also be found in the transfer- 
operator spectra discussed in Refs. @,[23|. Looking for an 
explanation for this similarity, we recall that the trans- 
fer operator is infinite-dimensional and independent of L. 
Thus, one may be tempted to compare Eq. with the 
infinite-L limit of the transfer matrix (|lip . In this limit, 
the transfer matrix becomes an integral kernel, indices 
becoming continuous variables: 



< g < 1 . 



(21) 



2'' - 1 

A careful calculation verifies that the transfer matrix 



5 




FIG. 4: (color online) Semiclassical eigenvalues (moduli) ver- 
sus scaled eigenvalue index k/N (decreasing order). Note 
the clear transition around k/N = 1. Parameters are A*' = 
82 and L = 9, 10, 11, 12 (black/solid, red/dash, green/dot, 
magenta/dash-dot, respectively). 



tends to the transfer operator, i.e., 

lim Mkk' = W{q, q') 



(22) 



with W{q,q') precisely that given in Eq. ([8])! So, the 
transfer operator is formally recovered as the continuum 
limit of the transfer matrix. 



IV. TRUNCATION SCHEMES 

Here we discuss the numerical calculation of long-time 
Gutzwiller traces using the transfer- matrix formula (|10p . 
The most natural procedure for calculating trAf^ re- 
quires the diagonalization of M to obtain its eigenvalues. 
This scheme is limited to very small values of L, e.g., to 
L ^ 15 in a personal computer. So, truncation becomes 
essential. 

The error introduced in truncating the interactions to 
r neighbors can be roughly estimated from Eq. ^ and 
the basic definitions (|5l6l7p . Such a truncation produces 
an error in the actions Q approximately given by 



SS ^ - e — e . 

2 2'- 



(23) 



For typical vectors e, containing randomly distributed 
0/1 bits, we have 



Finally, if we assume that this is the error committed in 
most phases in the Gutzwiller sum (O, then d(j) amounts 
to the relative error in the semiclassical traces, for 



(26) 



Thus, we arrive at the following condition for the validity 
of the truncation to r neighbors: 



(27) 



However, even after truncating the transfer matrix, we 
may need to use values of r which make diagonalization 
prohibitive, e.g., for N = 50, L = 1000. Remarkably 
there is an alternative to diagonalization which arises 
from the sparsity of M. For, even if the special type of 
sparsity M exhibits is not helpful in speeding-up its di- 
agonalization, it permits us to implement an alternative, 
much more efficient method for calculating the traces of 
M. 

Consider the following identity: 



trAfi 



2'- (V^IA/Z/^IV) 



(28) 



where overline means average over random complex vec- 
tors \ip) uniformly distributed over the corresponding 
unit sphere [3l|. (The notation Mr has been used to 
indicate the truncation of Af to r qubits.) The idea be- 
hind (pS)) is to calculate (tplMrltp) by applying iteratively 
the matrix Mr to \tp), and then averaging over This 
method allows one, in principle, to consider r as large as 
20. 

Table U displays some examples. The first rows (a-e) 
illustrate the improvement of the results as the trunca- 
tion level r is increased, for fixed {N, L). Agreement with 
the exact results is met, within the specified statistical 
errors, for r > 12, corresponding to / ^ 4. In rows (f- 
i) we consider (N, r) fixed and large values of L which 
can in no way be reached by direct computation of the 
periodic orbit summation, so that exact results are not 
known for such L. In these cases we compare with the 
approximate results obtained using the transfer-operator 
method 

We see good agreement for L < 80, i.e., we find the 
condition / ^ 4 again. 

So, we have checked that the transfer-matrix method 
works satisfactorily in the parameter regime specified by 

(123. 



DIVERGENCE OF LONG-TIME TRACES 



This implies that the corresponding error in the 
Gutzwiller phases amounts to 



S {2ttNS) 



NL 



(24) For chaotic systems, the Gutzwiller trace formula is 

only the first term of an expansion in powers of h. Be- 
cause of this, the semiclassical eigenvalues obtained from 
the Gutwiller traces deviate in general from the unit cir- 
cle in the complex plane. For typical chaotic systems it 
is expected that the distance of the largest semiclassi- 

(25) 

cal eigenvalue to the unit circle scales like fi, i.e., 1/N 
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L 


r 


TMM 


TOM 


exact 


a 


20 


10 


3. 73 [5] 


3.87[6] 


3.85 


b 


20 


11 


3. 90 [7] 


3.87[6] 


3.85 


c 


20 


12 


3. 85 [3] 


3.87[6] 


3.85 


d 


20 


13 


3. 84 [5] 


3.87[6] 


3.85 


e 


20 


14 


3. 89 [5] 


3.87[6] 


3.85 


f 


40 


13 


2.04[2]10^ 


1.99[1]10^ 


? 


g 


60 


13 


1.00[5]10^ 


1.05[5]10^ 


? 


h 


80 


13 


1.30[2]102 


1.24[1]10^ 


? 


i 


160 


13 


LOOfllW 


9.0[3]10=' 


? 



TABLE I: Gutzwiller's traces (absolute value). We com- 
pare the transfer-matrix method (TMM) with the transfer- 
operator method (TOM) and exact results. In all cases the 
quantum dimension is A'^ = 24. The square brackets con- 
tain the estimated error in the least significative figure [e.g., 
the notation 1.30[2]10^ stands for (1.30 ± 0.02) x 10^]. The 
values in the TMM column are the result of averaging over 
more than 10'' random states; the errors are statistical. In 
the TOM case the errors correspond to the variation of the 
results when the size of the Fourier basis is varied from 149 
to 201 [i,[23|. 

for quantum maps fs^]. In the case of the baker map, 
an anomalous behavior N~^/'^ is observed which is due 
to diffraction effects originating in the discontinuities of 
the mapping [IJ] . In any case, this lack of unitarity 
makes Gutzwiller's traces exponentially growing with L, 
a behavior that may be guessed from Table HI 

The following question arises naturally: can the N^^^^ 
law be associated with the particular structure of the 
baker transfer matrix? In the circuit representation 
(Fig. [2]) nonunitarity arises from the single-qubit gate H. 
So, one may ask: what is typically the largest eigenvalue 
(in modulus) of a map obtained by tensoring H with a 
generic A^-dimensional unitary gate U? 

On the other side, the g-representation ([TT]) offers a 
different point of view of the nonunitarity: the transfer 
matrix can be split as M = {U1+U2)/ V^, with Ui unitary 
matrices of dimension N [see Eq. ([31]) ]. How does the 
leading eigenvalue of such a matrix M scale with A^? 

In order to answer the questions above we resorted to 
numerical calculations. First we analyzed the matrices 
H' (El U with U a random matrix belonging to the Cir- 
cular Unitary Ensemble (CUE). The matrix H' acts on 
the least significant qubit and is defined by Eq. (IT71) with 
X = 1, i.e., H' is chosen to be the large-L limit of H. 
Then we considered the ensemble {Ui + \J-2)l\f2 with Ui 
independent random matrices belonging to CUE 33] . In 
order to check any possible ensemble dependence we also 
calculated the case H' ® U with U in the Circular Or- 
thogonal Ensemble (COE). 

Our results are exhibited in Fig. [5] In the three con- 
sidered cases we observe very similar decay laws. These 
decays are much closer to the diffractive scaling N~^^'^ 




FIG. 5: Average distance of the largest eigenvalue to the unit 
circle versus dimension for three ensembles of nonunitary ran- 
dom matrices. Triangles, squares, and circles correspond to 
H'®COE, (CUE-I^CUE)/^ (normalized sum of two indepen- 
dent CUE matrices), and iif'®CUE, respectively. See text for 
the definition of the single-qubit nonunitary gate H' . Error 
bars correspond to statistical errors. Straight lines represent- 
ing decay laws N~^^^, N~^^'^, and have been drawn for 
reference. 



than to the universal semiclassical behavior N~^, mean- 
ing that the random-matrix modeling has partially cap- 
tured the essence of the semiclassical baker. However, 
there is still some noticeable departure from the iV~^/^ 
decay (rather, our numerical results seem to follow a 
law). Thus, we conclude that the minimum- 
information models we have constructed still have to be 
complemented with some ingredient to properly describe 
the N^^/"^ baker scaling. Further research is necessary to 
discover what such an additional information should be 
(e.g., some particular correlation between Ui and i72?). 



VI. REFLECTED BAKER MAP 

The quantum baker map B is invariant under the par- 
ity symmetry i?, its action on the q-basis being just 
R\k) = \N — k). This is the quantum counterpart of the 
classical reflection symmetry TZ : {q,p) — {1 — q,l — p). 
When analyzing spectral properties it is convenient to 
separate eigenstates and eigenvalues of B according to 
their parity. Thus, one is led to consider the parity- 
projected bakers B± = {B±BR)/2. In the semiclassical 
domain the traces of B± can be approximated by Q 

1 / 9-f'/2 9^/2 , \ 

± 2 \^2i - 1 ^ 2^ + 1^ J 

(29) 

The first summation above correspond to the baker traces 
dS]). In the second sum 5"(e) stands for half the action 
[Eq. ^] of a periodic orbit of length 2L with symbolic 
code e = (eq, ei, • ■ • , fi-i, eb, ei, • ■ • , el-i)*, where e"i = 
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1 — Si. These are precisely the orbits invariant under the 
parity transformation. 

We shall exhibit the transfer matrix associated to the 
reflected traces, pointing out the most important fea- 
tures. Following the same procedure as in Sec. we 
constructed a transfer matrix M' such that 



2L 



(30) 



After eliminating the null subspace that appears as a con- 
sequence of considering just the parity-invariant trajecto- 
ries of length 2L, the matrix M' is reduced to dimension 
2^ X 2^: 




M' 



1 ° 


1 





••• • 
















1 • 

















••• • 





1 


ao 








••• • 














fll 


••• • 















••• 






V 







••• • 


•• a2L-i_ 


-1 0/ 



(31) 



where the matrix elements are given by 

~k 



ak = exp 



1 



with 



TiN2 



22L-1 



2L-1 



22i - 1 



(32) 



(33) 



The last undefined ingredient is the integer k. It is best 
expressed in terms of the binary digits of the integer k 
[see Eq. (IT5|) ]. If fc = e2i.-i_i . . . eiCo, then 



k = 621 



, .eieoOe2£ 



. .eieo . 



(34) 



The main difference between the baker's transfer ma- 
trix ([TT|) and the matrix above, is that M' is exactly 
unitary (thus, its spectrum lies on the unit circle). Ex- 
cept for this fact, both M and M' are structurally very 
similar. Perhaps this can be better appreciated by look- 
ing at the quantum circuit for M' in Fig. [S] The circuit 
acts on L qubits (without truncation). It starts like the 
baker's, with a downwards shift of all qubits. Then we 
have a NOT gate acting on the first qubit given by 



N 




(35) 



A phase gate -Poo (7) follows that acts on the first qubit, 
with 



7 = a' 1 



1L-1 



- 1 



22L-1 



(36) 



FIG. 6: Quantum circuit associated with the semiclassical 
traces of the parity-reflected baker map. From left to right, 
the circuit is composed of a downwards shift of all qubits {S^), 
a NOT gate (A''), a one-qubit gate, and a sequence of L — 1 
two qubit phase gates. All gates are unitary. 



Finally there is a sequence of L — 1 symmetrical two- 
qubit phase gates Pofc(/^) acting between the first and 
the fc-th qubit with exponentially decreasing phases, j3 = 
5/2, S/4, . . . , 5/2^, where a = a' {2^ - l)/2^. 

A complete study of the spectral properties of the 
reflected-baker transfer matrix, truncation schemes, etc, 
is deferred to a future publication because, as M' is uni- 
tary, such analyses would take us in directions very dif- 
ferent than those followed in the case of the baker map. 
However, we would like to anticipate one result which 
gives a hint about the nature of the spectrum of M' . 

Consider the 2L-th power of the transfer matrix M', 
i.e., the 2L-th iteration of the circuit above. We verified 
that the shifts and NOT gates cancel out and only phase 
gates remain. This means that M''^^ is a diagonal ma- 
trix (given that the phase gates are diagonal). It turns 
out that the diagonal elements are exactly the Gutwiller 
phases exp{iS' /h) in Eq. pS)) (we tested this numerically 
for some cases). Thus, the spectrum of M' is formed by 
roots of the Gutzwiller phases! 



VII. CONCLUSIONS 

We presented a study of the transfer matrix approach 
to the semiclassical traces of the baker map and its re- 
flected version. We found that the transfer matrices ad- 
mit a tensor product decomposition leading to simple cir- 
cuit representations, similar to those obtained by Schack 
and Caves for the quantum baker. Remarkably, in the 
case of the baker, the corresponding circuit clearly iso- 
lates the source of nonunitarity of the semiclassical traces 
in the form of a one-qubit Hadamard-like gate. Given 
that both exact and semiclassical bakers can now be writ- 
ten as circuits, this representation appears as a promising 
tool for studying the corrections to the Gutzwiller trace 
formula. 
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In the case of the baker, spectral properties were an- 
alyzed, showing that there exists a close similitude with 
the transfer operator of Ref. @. In fact, we proved that 
the transfer operator arises as a suitable asymptotic limit 
of the transfer matrix formalism. Truncation schemes 
were discussed which permit the numerical calculation of 
long-time traces well inside a domain where direct sum- 
mation of the Gutzwiller formula is impossible. 

We would like to conclude by mentioning a very excit- 
ing (though rather speculative) possibility. It is tempt- 
ing to think of the transfer matrices as maps for which 
the Gutzwiller formula is exact. Admittedly this point 
of view has some limitations because there is a map for 
each value of L. Even so, our results suggest that we 
may be not far from finding a close relative to the baker 
map having an exact periodic-orbit trace formula. The 



search for such a map constitutes a very attractive and 
challenging project which, however, exceeds the scope of 
the present paper. 
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